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Abstract 

Here we estimate the magnitude of the signals expected for realistic 
cosmic neutrino backgrounds in detectors attempting to measure the me- 
_C ' chanical forces exerted on macroscopic targets by the elastic scattering of 

CIh, relic neutrinos. We study effects proportional to the weak coupling con- 

stant Gf and to G F for Dirac and Majorana neutrinos, either relativistic 
qj ' or non-relativistic, both gravitationally bound or not. 



1 Introduction 

^h ' Alongside the well measured background of cosmic photons, the Big Bang model 

predicts the existence of an elusive background of cosmic neutrinos. This neu- 
trino background may be revealed in the near future by measurements of the 
microwave background radiation anisotropy and the large-scale structure of the 
Universe. 

Here we reconsider the long-standing question of the detectability of the 
neutrino background in a laboratory experiment, a recurrent subject in the lit- 
erature for the past thirty years fl|-[f§. We estimate how far macroscopic accel- 
erations due to realistic neutrino backgrounds are from the smallest measurable 
acceleration at present, which is about 10~ 12 cm/sec 2 0. The new elements we 
consider in our estimates are the possibility of a very large lepton asymmetry 
in the cosmic neutrino background, and the study of Majorana as well as Dirac 
neutrinos. 



With no lepton number asymmetry the number density n Vi of the neutrinos 
of each species in the background is equal to the number of anti-neutrinos of 
each species. The Big Bang model predicts this number density to be n v = 
(3/22) n 7 = 0.136 n 7 , where n 7 is the density of the microwave background 
radiation photons. Since n 7 = 412 cm -3 , therefore n Vi = 56 cm" 3 . It is, 
however, possible for neutrinos to have very large lepton asymmetries. While 
charge neutrality requires the asymmetry in charged leptons to be the same 
as that in protons, for which (ng — n^)/n 7 ~ O(10 -10 ), no such requirement 
limits the asymmetry in neutrinos. With a large asymmetry between neutrinos 
and anti-neutrinos only the excess of either Vi or Vi remains after annihilation 
ceases with a density 

n Vi = n i 0.0252 [&7T 2 + ef] (1) 

where £j = \/-i Vi \/T Vi and \x Vi is the chemical potential of the given neutrino 
species. The ratio £j is constant after neutrinos decouple. 

Thus, for example, with £, = 5 [pi O, QQ], one obtains a neutrino density of 
n v ~ 4n 7 = 1700 cm -3 . Eq. (1) is valid for ^ < 12 || (for which the decou- 
pling temperature Tdec is smaller than 2771^, see below). Bounds coming from 
Big Bang Nucleosynthesis (BBN), the anisotropy of the microwave background 
radiation and data on the large scale structure of the Universe limit £j to be 
much lower than 12. Kang and Steigman in 1992 || obtained 

- 0.06 <£ e < 1.1 , |^,r|<6.9. (2) 

More recent studies yield bounds that depend on the presence of a cosmological 
constant (with energy density Q& in units of the critical density) . If the cosmo- 
logical constant is zero or small, £7 a < 0.1, the value £j ~ 5, namely n Vi jn 1 ~ 4 
is favored P , jlCfl ■ For CIa — 0.5, one obtains n Vi /n^ < 4, and for Q\ > 0.7, 
the lepton asymmetry must be very small, so the neutrino density goes back 
to the standard value, n Ui /n-y ~ 0.14 ||, ^l). Notice that in models with large 
neutrino asymmetries, the standard BBN mechanism is altered, allowing for the 
existence of many more neutrinos than in the standard case. BBN and recent 
cosmic microwave background anisotropy data provide constraints on neutrino 
asymmetries which depend strongly on cosmological parameters JL2| . In what 
follows we will take an upper bound n Vi /n 1 < 4, i.e. n Ui < 1700 cm -3 , when 
considering the case of large lepton asymmetries. 

The decoupling temperature of neutrinos T^cc, which with no lepton asym- 
metry is a few MeV, increases with the lepton asymmetry. However, only for 
li Vi lT Vi > 12 would Tdec become larger than 2m ^, i.e. neutrinos would decou- 
ple before fi + [i~~ annihilation. In this case, not only would e + e~~ annihilations 
increase the number of photons relative to that of neutrinos due to entropy 
conservation, but /i + /i~ annihilations would as well, leading to a lower neu- 
trino temperature relative to the photon temperature. For fj, Ui /T Ui < 12 the 
temperature of Pi is now the usual one, T Vi — (4/11) 1 / 3 T 7 = 0.71 T 7 . With 
T 7 = 2.728°K, this means T Vi = 1.95°K =*1.68 10~ 4 eV. Therefore, if m Vi is 



smaller than 1.7 10 -4 eV background neutrinos are relativistic at present while 
for larger masses they are non-relativistic. Moreover neutrinos may be Dirac 
or Majorana particles, a distinction that is important only for non-relativistic 
neutrinos. 

Relativistic neutrinos are only in left-handed chirality states (and anti- neu- 
trinos only in right-handed chirality states). These are the only states produced 
by weak interactions. For relativistic neutrinos chirality and helicity coincide 
(up to mixing terms of order m Vi jE v ~ ra Vi jT v ). Bounds from structure forma- 
tion in the Universe imply that stable neutrino masses are at most of the order 
of a few eV [|3[. This means that neutrinos were relativistic at decoupling 
(?dec > O(MeV)), even if they may be non-relativistic at present. 

We will call left(right) chirality eigenstates vl{vr) and left(right) helicity 
eigenstates vi{v r ). At decoupling, neutrinos ul were only in left-handed helic- 
ity states and antineutrinos i/ R (or vr in the case of Majorana neutrinos) in 
right-handed ones. Helicity is an eigenstate of propagation and therefore it does 
not change while neutrinos propagate freely, even if they become non-relativistic. 
For Majorana neutrinos chirality acts as lepton number, so we are calling "neu- 
trinos" those particles produced at T > T^ ec as vl and "anti-neutrinos" those 
produced as vr. Thus, neglecting intervening interactions, non-relativistic back- 
ground neutrinos are in left-handed helicity eigenstates (which consist of equal 
admixtures of left- and right-handed chiralitics) and anti-neutrinos are in right- 
handed helicity eigenstates (which also consist of equal admixtures of left- and 
right-handed chiralitics). If the non-relativistic neutrinos are Dirac particles, 
only the left-handed chirality states (right for anti-neutrinos) interact, since 
the other chirality state is sterile, while if the neutrinos are Majorana, both 
chirality states interact (the right-handed "neutrino" state is the right-handed 
anti-neutrino). 

Slow enough non-relativistic neutrinos eventually fall into gravitational po- 
tential wells, become bound and, after a characteristic orbital time, their he- 
licities become well mixed up, since momenta are reversed and spins are not. 
Thus, gravitationally bound background neutrinos have well mixed hclicitics. 
Most background neutrinos however are not gravitationally bound at present 
because they are too light. 

The present upper bound from structure formation in the Universe for neu- 
trinos without a lepton asymmetry, for which n„ ~ 100 cm~ 3 , is of a few 
eV [O. Neutrinos carrying a large lepton asymmetry must necessarily be 
lighter. The upper bound on the mass of degenerate neutrinos with a large 
asymmetry should be lower, since their number density is larger, at most of 
the order of 0.1 eV. The Tremaine-Gunn kinematical constraint |lj] gives the 
scale at which neutrinos gravitationally cluster. Light neutrinos with masses 
m v < eV would be gravitationally bound only to the largest structures, large 
clusters of galaxies. We can see this using simple velocity arguments. Only 
cosmic neutrinos with velocities smaller than the escape velocity of a given 
structure can be bound to it. The escape velocity from a large galaxy like 



ours is about 600 km/s and from a large cluster of galaxies is about 2,000 
km/s. Considering that the average velocity modulus of non-relativistic neutri- 
nos of mass m and temperature T v is (using Maxwell-Bolztman distribution) 
(\fj u \) = v / 8kT u /nm= v / 4:.3 W^eV/m (namely (\v v \)= 6,200 km/s for m = leV, 
and (|^|)= 19,600 km/s for m = O.leV) it is obvious that only about a third 
of 1 eV mass neutrinos and a very small fraction of lighter neutrinos could be 
gravitationaly bound to large clusters at present. Fermi degenerate neutrinos 
may have even larger average velocities depending on their chemical potential 
£ = fi/T v (constant after neutrinos decouple), but the conclusions remain the 
same. For £ » 1 (\J3 v \) = y/6£T v /5m~ a/£1.68 lO^eV/m (namely (\v v \)= </£ 
12,300 km/s for m = O.leV) and both expressions coincide for £ = 2.5. In all 
cases the amount of neutrinos in the tail of the velocity distribution with veloc- 
ities smaller that 600km/s which would be gravitationally bound to galaxies is 
much smaller. 

In the following we will give our results for both clustered (C-NR) and non 
clustered (NC-NR) non-relativistic background neutrinos as well as for relativis- 
tic (R) background neutrinos. 

With zero lepton asymmetry the background consists of equal numbers of 
neutrinos and anti-neutrinos. Thus we have equal numbers of left and right 
either chirality or helicity states in all cases, relativistic or non-relativistic, Dirac 
or Majorana neutrinos. In the case of a large lepton asymmetry which favors, 
say, neutrinos so that there are only vl at the moment neutrinos decouple, 
(left and right are to be exchanged in the following argument for an asymmetry 
favoring anti-neutrinos) the background consists of left chirality particles in the 
case of R or NC-NR neutrinos, and it consists of totally mixed helicity states 
for C-NR neutrinos 

These properties of the different possible neutrino backgrounds are relevant 
to the effect linear in the weak coupling constant Gp, first studied by Stodolsky 

0- 

2 The G F Effect 

This effect was proposed by Stodolsky in 1974 (before weak neutral currents were 
proven to exist, so he did not include them). It is the energy split of the two 
spin states of non-relativistic electrons in the cosmic neutrino bath. This energy 
split is proportional to the difference between the densities of neutrinos and 
antineutrinos in the neutrino background for Dirac neutrinos, and proportional 
to the net helicity of the background for Majorana neutrinos, as we will now 
see. 

The Hamiltonian density of the v-e interaction is 

C 
U{x) = — = ej^igv - 9Als)e v lti (\ - l5 )v . (3) 



For v = v e , qa = 1/2 and gv = 1/2 + 2 sin 9w while for v = v^ or v = v T , 
with which electrons only have neutral weak interactions, gA = — 1/2 and gv = 
-l/2 + 2sin 2 6> w . 

2.1 Dirac Neutrinos 

We compute first the energy shift of a single electron of momentum p and spin 
s in the neutrino background: AE e — (p,s\H\p,s). Let us call p and k the 
momenta of the e and v respectively, and s(s') the incoming(outgoing) spins. 
Working in momentum space and using the second quantization form for the 
neutrino fields, we obtain for the case of Dirac neutrinos 
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where it and u are the usual spinors for particles and antiparticles and 6, d 
(b\ d') are the respective annihilation (creation) operators, for neutrinos and 
antineutrinos respectively, of momentum k. In the rest frame of the electron 
(p = 0) the diagonal matrix element of the weak current is gA^l- In this frame 
the neutrino current is 3 % v — — (3 l earth (n v — n v c), since the non-zero average 
velocity —j3 ear th of the neutrinos is due to the motion of the earth relative 
to the neutrino bath (with velocity dearth)- Dotting the two together yields 
AE = -(G F /\/2)g A (a e ■ (3 eart h) {n v -n v c). 

We would like to present however a more careful, systematic, and first 
principle derivation. Using the decomposition u^(p, s) = \/(E e + m e )/2m e 
(Xs '> [(^e • P)/(E e + r7 ie)]xJ 1 ) f° r the electron spinors and the Gordon decom- 
position and similar relations containing 75 for the neutrino vertices, one obtains 
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Here the comma in the first key bracket separates the zero and spatial com- 
ponents of the electron amplitude and W^ = i<J^ v k v ^ = —^e filJap cr a ' p k L ' is the 
Pauli-Lubansky pseudo- vector (one uses here 750^1, = |e /iyCTp tT crp ). In terms of 
the four dimensional "spin operator" 



£ = 



a 
a 



the Pauli-Lubansky pseudo-vector is given by 



E, 



W »= o - fc > -^fS + - 75 (/cxS) , 



(6) 



(7) 



so that for particles at rest (fc = 0) it is proportional to the spin operator 
Wa = (0, — m„S/2) and for relativistic particles it is proportional to the 
helicity h = S • k, i.e. W^\k) — (h/2)k ll \k). Although W is the spin operator in 
the rest frame and W^W^ — —m 2 s(s + 1) is a Casimir operator of the Poincare 
algebra, in general W is not the spin operator: it is easy to see that the Wi 
do not have the SU(2) commutation relations required for angular momentum 
components, i.e. [Wi, WA 7^ ieijkWt- However, one can write W in terms of 
the true spin operator |l5[ 



m v \ E u + m v 



(8) 



Considering non-relativistic electrons with velocity [3 e = p/E e (thus (E e + 
m e )/E e — 2), using (3 V — k/E v , and writing W in terms of S, we obtain to 
first order in f3 e 

(H D ) = -j= {(gv ~ 9A 2s e ■ e ), 2s e (-g A + g v 2s e • /3 e ) J 
7 I fd 3 kJ2b\k,s')b{k,s) \u v {k,s')(l-t-p v , 

I. S i S> 
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-fl + ^S+ J^" g • & ] v„(k,s') 



(9) 



In the neutrino "rest-frame" the term of the form s e ■ S in Eq. ffl) vanishes 
for the following simple reason. We have calculated our interaction between 
a "sea" of neutrinos and a single electron. Thus s e ■ S gives the projection 
of the individual neutrino spin along the single electron's spin axis. But we 
expect the neutrinos to be distributed isotropically, and thus there is equal 
probability for the neutrino spin projected along this axis to be positive or 
negative. Hence when the sum over spins and the integration over the isotropic 
neutrino distribution are performed, these terms becomes zero. However, in the 
electron rest-frame the neutrino distribution is no longer isotropic and this term 
does not vanish. Using the basis of neutrinos of definite helicity, we can take 
the spin to lie completelty parallel or anti-parallel to the direction of motion, 
i.e. S — (h/2) (3 where h is the helicity. Thus the term in Eq. (g) containing 
the neutrino spin will be proportional to ((m^/E^) f3). 

Instead of utilizing the spin basis, we now work in the helicity basis so that 
all spin dependent terms involve the helicity operator £ • /3. We can use now 
the completeness relations for the spinors u and v, 

u u (k,h')u„(k,h) = Shh' v u (k,h')v u (k,h) = -<W- (10) 

Furthermore, in the helicity basis we have 

u v (k, h!) Ej/ • (3 V u„(k, h) = h\(3 v \ Shh', (11) 

v v {k > h')t v -0 v v v (k,h) = -h\^\5 hhf , (12) 

where h = 2s-/3 is the helicity. Thus we obtain 
(H D ) =*~{( gv -g A 2s e ■ e ) , 2s e (-g A + g v 2s e ■ /3 e )} X 
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The number density of neutrinos and antineutrinos are respectively 

n v = — f d 3 kJ2 b Hk,h)b(k,h), (14) 

■* h 

- I d 3 k^2d^{k,h)d{k,h), (15) 

j i. 



and obviously 



V 



i J d*k^b\k,h)b{k,h)p v = 0u)n v , (16) 

h 



and similarly for antincutrinos, where ( ) denotes average values. 

Recall that we denote left- and right-handed helicity states with lower case 
indices, vi and v r . In general the neutrino and the anti- neutrino are admixtures 
of states of left and right helicity. Thus we write 

n v = n vt + n Vr , n v c = n v o + n K . (17) 

With this notation the terms in {H D ) linear in s e are H®, 

H Z = (G F /V2)l~g A 2s e - p e (n v -n u c) + (\/3 v \) (n Vl +n v c - n Vr - n„<=) 



+ 9A 2s e 



(Pv){n v -n u c)+ 

-Pv\Pv\ ) + ( -£r Pi> ) ) {n Ve + n v o - n Vr - n K ) 



E v +m v / \E V 



(18) 



In the rest frame of the i/-bath due to isotropy the average of all vectors pro- 
portional to neutrino velocities (/?„) = 0, (J3 U ) = 0, {{E v /E v + m u )(3 u \fl u \) = 
and ({m v /Ev) J3 y ) = 0, thus only the 1st. term in Eq. (18) remains. In the rest 
frame of the electron, the frame in which experiments are performed, j3 e = and 
only the 2nd. term in Eq. (18) remains. In this frame the non-zero value of (/?„} 
is due to the motion of the earth relative to the neutrino bath, (/3 W ) = —dearth, 
which is approximately 10~ 3 . One can compute the averages in Eq. (|18|) ex- 
plicitely in the relativistic (R) and non-relativistic (NR) limits, using a Gaussian 
distribution for the neutrinos in the latter case. Up to first order in (3 ea ,rth we 
have 

L " -/i|/l|\ =-{\jL\)0earth = -dearth, (19) 

/ R 

m^\ = _&o^ (2Q) 

Ev I R 7 
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^Pu) =-^^ + 0(Pl arth ), (22) 

E v / NR 3tt (I/J^l) 

where (\P U \) is always the average of the velocity modulus in the neutrino rest 
frame. We see in Eqs. (19) and (21) that {{E v /E v + m v )p v \p v \) is either of 
the same order of magnitude or smaller than p ear th- The remaining average, 
{(m v /E v ) p v ) is negligible for relativistic neutrinos with 7 >> 1, but is larger 
than Pearth for non-relativistic neutrinos, for which (|/3„|) < 1 (see Eq. (22)). 
The factor (n Vt + n v c — n Vr — n u c ) , which multiplies both these averages in Eq. 
(18), becomes zero for C-NR neutrinos (since the helicities are well mixed) but 



could be large for NC-NR neutrinos in the case of a large lepton asymmetry. 
In this case -{{m v /E„) 0„) = 1.7^/^1.7 lG~ A eV f3 earth < (14/V?)/3 e0 rtfc, 
(considering that the mass of neutrinos with a very large lepton asymmetry is 
at most of the order of 0.1 eV) and this term is dominant. 

^From Eqs. (18) to (22) we have that in the rest frame of the electron, the 
frame in which experiments are performed, to first order in (3 ear th 



+ 



(x{\J3„\)- 1 + y{\0 v \)) (n vi + n v c - n Vr - n„.)] (23) 



where y=l and x = for relativistic neutrinos and y=2/3 and x = 16/3-7T = 1.7 
for non-relativistic neutrinos (these numbers are the prefactors in Eqs. (|l9|), 
(pTj) , and ( |22[ ) above). While in the rest frame of the neutrino bath, where 
(P v ) =0 and (i e = f3 earth we find 

C 1 - r - i 

Hf e = — j= 9a 2 s e ■ p e arth (n v - n v a ) + ( \ j3 v \ ) {n vi + n„o - n Vr - n v o ) . (24) 

For relativistic neutrinos, for which chirality and helicity coincide, we find 
in both reference frames that H® = H® is 

Hft = -V2G F 9a 2s e • [i e arth{n UL - n„c) = -V2G F gA 2s e • Pearth{n v - n u <=) , 

(25) 
and for non-relativistic gravitationally bound (C-NR) Dirac neutrinos (with well 
mixed helicities) we find to order (3 ear th, again in both reference frames, that 

H Z = H C-NB. is 

C 1 

H C-NR = — JE9A 2s e ■ Pearth{n v - n v c) = -Hr , (26) 



which is smaller than the effect in Eq. (|25j) for relativistic neutrinos by a factor 
of two. 

For most of the non-relativistic relic neutrinos, those non-clustered, in the 
presence lepton asymmetry, we find that in the rest frame of the electron the 
dominant contribution to H® = H® C _ NR is 

H NC~NR — -7^9A^S e ■ P<. a rthl-7{\Pv\)~ l {n Vl + Tl v c - n Vv - 7l„c) , (27) 



with (IA,!)- 1 = y/m/Z1.7 10- 4 el/ (3 earth < (8.2/^){3 earth . 

It is obvious that the effect is non-zero only in the presence of a lepton 
asymmetry, where n u ^ n v c . The effect is maximum if the relic bath consists 
only of neutrinos (or only of antineutrinos) so that n v ^ and ro„c = (or 



viceversa), which is possible with a large lepton asymmetry (fi/T > 2). In this 
case, Eq. (p5|) becomes 



H% = -V2G F g A 2s e ■ (3 earth n u , (28) 

and Eq. (j27|) becomes 



Hgc-s* - - 8 VotTFw^ * ^ < (29) 

since m v < O.ley, in the presence of a large lepton asymmetry. 

2.2 Majorana Neutrinos 

In comparison to Dirac neutrinos, Majorana neutrinos satisfy an additional 
constraint, v = v c . We write the Majorana field as ** f = (^ + (*^) C )/V2. 
Therefore, using the ordinary decomposition for Dirac neutrinos *§>„ one arrives 

at 
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(30) 



Defining a new operator 6(fc, s) = (6(fc, s) + d(k, s))/v2 the general wave ex- 
pansion for a Majorana field is given in terms of only one creation and one 
annihilation operator, b and w . Note that the factor l/\/2 is included in the 
definition of b and w , instead of keeping it as an overall factor. The overall 
normalization of a Majorana field is trickier than that of a Dirac field. Our 
Dirac fields are normalized such that J d 3 x^^ ty® — N y — N„. However, in 

the case of Majorana fields, J d 3 X^fif *$!„ — 0, and obviously such a condition 
cannot be used to determine an overall normalization. However, one can check 
the anti-commutation relations of the creation and annihilation operators b and 

w . One finds that < b(k, s), w (k, s') > = 5 s , S 'S 3 (k— k'), which is the correct rela- 
tion and proves that the definitions of 6 and b' are correct and that no further 
normalization of ^* f is needed. 

Since the general wave expansion of the Majorana field is written in terms of 
only one creation and one annihilation operators, one can recover the result for 
Majorana neutrinos directly from the calculation for Dirac neutrinos by taking 
b = d in Eq. (^|). Effectively the vertex 7^(1 — 75) in Eq. (^) becomes a pure 
axial 7^75 vertex. This amounts to the absence of the fc M terms in Eq. (g) and 
thus the terms proportional to (n v — n„<=) in Eq. ( |l8| ) are absent. Therefore, 



10 



the effect depends on having a net helicity in the Majorana-neutrino bath. The 
terms linear in s e in H are now: 



H™ 



V2G F g A {-2s e -p e {0„\) 



In the rest frame of the electron, the frame relevant for experiments, this is 

H^ = -V2G F g A 2s e -p earth [x{\p v \)- l +y{\f3 v \)\ (n vi - n Vr ) (32) 

where, as before, x = and y=l for relativistic and a;=16/37r and j/=2/3 for 
non-relativistic neutrinos. For relativistic Majorana neutrinos, for which helicity 
coincides with chirality vt = vl and v r = Vj>, this term coincides with the result 
for relativistic Dirac neutrinos in Eq. (Eq), as it should be since Dirac and 
Majorana neutrinos cannot be distinguished when relativistic. In the absence 
of interactions the helicity is conserved, even when as the temperature of the 
universe decreases neutrinos become non-relativistic. In the case of a large 
lepton asymmetry favoring, say, neutrinos, so that only v^ = v^ are present at 
the moment of decoupling (when neutrinos are relativistic) a net left helicity 
remains in the bath of Majorana neutrinos, while not gravitationally bound. In 
this case, without gravitational binding, Majorana neutrinos would have a net 
left helicity n vi = n VL 7^ 0, n Vr = n u c = 0. As in the case of Dirac neutrinos 

the term proportional to (|/3i/|) _1 dominates for non relativistic neutrinos. This 
is important only for unclustered (NC-NR) neutrinos. As already mentioned, 
gravitationally bound non-relativistic (C-NR) neutrinos loose their net helicity. 
After a characteristic orbital time, helicities become maximally mixed thus I 
and r helicities become equally abundant, so (n Vl —n Vr ) — 0, and the Stodolsky 
effect vanishes. 

Thus, to summarize, in general 

K=K (33) 

^C-NR = (34) 

and for neutrinos with a large lepton asymmetry, so that say n Ul , ~ n„, 



In this last equation H^ is that given in Eq. (|2q). 

2.3 Accelerations 

The Stodolsky effect consist of a difference in energy between the two helicity 
states of the electron interacting with the neutrino bath in the rest frame of 



11 



the earth in which experiments are performed (which we take to coincide with 
the rest frame of the electrons) . This energy difference AE is obtained in each 
case by replacing s e ■ f3 ear th by \/3 e arth\ in the respective expressions for the 
H D and H M given in the previous sections. Therefore, the Stodolsky effect 
vanishes except in the presence of a lepton asymmetry which in the case of 
Majorana neutrinos persists in the neutrino bath as an asymmetry in hclicity 
as long as neutrinos are not gravitationally bound. In the case of relativistic, 
light neutrinos of density n„, with a very large lepton asymmetry favoring, say, 
neutrinos v^, so that n UL = n Vl = n v , from Eqs. (pq) and (p3|) one has 



(A£)g = (AE)% = 2V2G F g A \/3 earth \n„ . (36) 



In the case of non-relativistic, unclustered neutrinos, we see from Eqs. (29) and 
© that 



(AE)% C _ NR = 2{AE) D NC _ NR ~ l-7 ^ 17 7 -4 e y ( A ^)g ^ ^ AE ^ ■ 

(37) 
In the case of gravitationally bound (C-NR) neutrinos (so that there is no net 
hclicity in the bath) of density n vi if there is a very large lepton asymmetry 
favoring neutrinos Vl, for Dirac neutrinos n„c = 0, and for Majorana neutrinos 
n Vl = n Vr and the effect vanishes. Thus, 

(AE)%_ NR ~ V2G F g A \J3e a rth\n u , (38) 

(AE)%_ NR ~ . (39) 

Equivalent results would be obtained with an asymmetry favoring antineutrinos. 

^From now on, we will use the energy difference AE in Eq. ([36|) to compute 
the maximal possible strength of the Stodolski effect, recalling that the effect 
could be at most one order of magnitude larger for non-clustered non-relativistic 
(NC-NR) neutrinos (most of the relic neutrinos if they are non-relativistic). 

The difference in energy AE between the two helicity states of an electron in 
the direction of the bulk velocity of the neutrino background < f3 v >= —f3 ear th 
implies a torque of magnitude AE/tt applied on the spin of the electron. Since 
the spin is "frozen" in a magnetized macroscopic piece of material with N 
polarized electrons, the total torque applied to the piece has a magnitude 
r = NAE/n. Given a linear dimension R and mass M of the macroscopic 
object, its moment of inertia is parametrized as I = MR /j, where 7 is a 
geometrical factor. In the typical case of one polarized electron per atom in 
a material of atomic number A, the number N above is N = (M/gr )Nav/A 
(using cgs units), where Nav is Avogadro's number. Thus, the effect we are con- 
sidering would produce an angular acceleration of order a = t/I and a linear 
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accelerations of order a,G F = Ret in the magnet, given by 

^av Aff 7 

where the Gp subindex indicates the mechanism we have described. 



Using the expression in Eq. (36) for AE with n v — / 100cm 3 we then find 

m-27ir CY\ f 100 \ f Cm \ ( Pearth\ Cm 

*Q r = 10 / • (-) (^— j (— ) (^ TFF j ^ (41) 

where / accounts for a possible local enhancement of the standard background 
neutrino number density. This acceleration could be at most one order of mag- 
nitude larger for non-relativistic (NC-NR). These accelerations are rather weak. 

As mentioned above the present upper bound from structure formation in the 
Universe for neutrinos without a lepton asymmetry, for which n v ~ 100 cm -3 , 
is of a few eV. The upper bound on the mass of degenerate neutrinos with a 
large asymmetry should be lower, since their number density is larger. Light 
neutrinos with masses m v < eV would be either unbound or gravitationally 
bound to very large structures and thus local enhancement to the neutrino 
density due to gravitational clustering would be very small. 

Throughout we have used /3 e arth ~ 10~ 3 , which is the velocity of the earth 
with respect to the galaxy; however, if neutrinos are gravitationally clustered 
on much larger scales, one must consider the relative motion of the earth with 
respect to the rest frame of these larger objects. At such scales the neutrino 
bath rest frame can be taken to coincide with the cosmic microwave background 
(CMB) rest frame. The sun's motion with respect to the CMB is believed to 
be responsible for the largest anisotropy in the COBE DMR maps, the 3 mK 
dipole, and thus has been determined with great accuracy. From the COBE 
data, v sun = 369.0 ± 2.5 km/sec which corresponds to (3 ear th — 1-231 x 10 -3 
[ pTG| . Therefore, even if background neutrinos are bound to such large objects 
as super-clusters, we are justified in using still (3 ear th ~ 10~ 3 . 

Thus with very little gravitationally clustering locally, a neutrino density 
enhancement could only be due to the asymmetry itself, in which case (as ex- 
plained in the introduction) / < 20 (n V( < 1700 cm -3 ). 

However, it is amusing to note that if there is a local cloud of v e of density 
~ 10 17 cm -3 as required for a rather audacious explanation of the tritium 
end point anomaly [|l7fl then aa F — 10 -12 cm/sec 2 , and the effect would be 
observable! 

2.4 Comparisons with Past Results 

The first calculation of this order Gf effect was done by Stodolsky in 1974 JJJ. 
In our notation, he found the energy difference between the two spin states of 
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an electron moving with velocity /3 e in the rest frame of the neutrino bath to be 
A£ = 2V2G F -jh (n v -n v c) . (42) 

V 1 - Pe 

Stodolsky performed his calculation before the discovery of neutral currents 
and thus did not include them (hence he took §a to be 1). He considered only 



relativistic neutrinos. Stodolsky's result conincides with ours in Eq. ( |36|) with 



f3e = \Pearth\ (and y/T^f 2 « 1). 

Langacker, Leveille, and Sheiman also explored this effect H. Although the 
majority of their paper deals with exposing the flaws in several relic neutrino 
detection schemes, they conclude that the only legitimate effect of order Gf is 
that proposed by Stodolsky. They studied only Dirac neutrinos and found 

AE = 2V2G F [3 e N t a ot ^ (n Vi - n„«) g%K(p, m t ) , (43) 

i 

where -/V*°* is the total number of electrons in the sample, and K is a function of 
E v and m v which reduces to 1 if m„ = or 1/2 if k v -C m v . Their calculation 
is performed in the rest-frame of the "neutrino sea" . Their result agrees exactly 
with Eqs. ( p6| ) and (j^) in which we give the AE for Dirac neutrinos in the rest 
frame of the electron. 

Ferreras and Wasserman g recently tackled the subject of the detectability 
of relic neutrinos, and also concluded that there are no order Gf effects without 
lepton asymmetry. However, they point out that density fluctuations in the the 
v(v c ) background could give rise to non-steady forces of order Gf- Such forces 
would give rise to displacements in massive objects proportional to t 3 ' 2 rather 
than t 2 as in constant acceleration. They caution that although such accelera- 
tions may be more readily detectable, they could also act as an additional source 
of noise. Here we have not considered neutrino density fluctuations. 

3 The G\ Effect 

Nucleons are continuously bombarded by the relic neutrinos and a momentum 
Ap « p v is imparted in each collision. The momentum of the neutrinos p u is 
p v = E v jc k, 4T u /c for relativistic neutrinos (R) and also for non-clustered, 
non-relativistic neutrinos (NC-NR) due to momentum red-shift, and is m v v v w 
muVyiriai for non-relativistic clustered neutrinos (C-NR). 

If the earth was at rest with respect to the relic neutrino "rest frame", 
i.e. the frame in which the neutrinos are isotropically distributed, then the 
average momentum transfer (Ap) would vanish. However the motion of the 
earth with velocity cj3 ea rth = v eart h induces a "dipole" distortion of 0(J5 ear th) 
in the velocity distribution of the relic neutrinos. In the laboratory frame (the 
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earth's rest frame) this makes 

(Ap) R « earth (E v /c) , 

( A P) NC -NR « dearth (4T„/ C ) = (Ap) fl , (44) 

( A P)c~NR ~ dearth C TO„ . 

The fluxes of infalling neutrinos are <£>_r = n v c and $jv.r = n v v v for rela- 
tivistic and non-relativistic neutrinos respectively. To find the resulting acceler- 
ations we compute the force exerted on one gram of detector material containing 
M =Nav /A nuclei (A,Z). The force is given by the momentum imparted per 
second: the latter is the microscopic Ap of Eq. (J4J) times NQ v <t v -a-, the 
number of collisions per second inside one gram: 

- = a = $ v ^a„- A {Ap). (45) 

m A 

The neutrino-nucleus cross sections <j v -a are extremely small even if we 
include a nuclear coherence enhancement factor (A — Z) S3 A 2 relative to the 
neutrino- nucleon ct„_jv cross sections, which are of the order 



_ ~ / G\<h *■ 10- 56 K {eV} ) 2 cm 2 for (NR) 

V N ~ \ G 2 F £> ~ 5 x 10- 63 (T. {1 . 9 o K} ) 2 cm 2 for (R) 

where { } indicates the units of the mentioned quantities. 
Collecting all terms in Eq. ( fl5| ) above we find 

l3 



Nau Gw 2 I ( 4T,y ) /3ear *' 1 for ( R ) 

'■ ' m„HT v p earth for(NC 

m v 3 v u 2 for (C-NR) 



//,. A 2 { m v 2 AT u (3 earth for(NC-NR) . (47) 



A " I 3„, 2 



We will take the earth's velocity v ear th and the clustered neutrino's velocity 
to be all ss cf3 V i r iai with f3 V i r iai ~ 10~ 3 corresponding to typical galactic virial 
velocities. We use n v = / 100 cm~ 3 . Hence we find 



KA JT- / \ 3 Ol 

a « = 3x10 /^(V«})^. ( 48 ) 



/I O CTT) 

0-nc-nr = 0.6 x 10~ 47 /— - (m„ {ev i) T^/i.gojf} g • ( 49 ) 

100 sec 2 

Here the density enhancement factor / may only be due to a large lepton 
asymmetry and / < 20 (see Sect. 1). Finally for clustered non-relativistic 
neutrinos we have 

ac-NR = io- 46 /^ Km) 3 ^ • (so) 

100 l '' sec 1 
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Here f contains also a clustering enhancement factor. Values of / up to 10 7 
have been mentioned in the literature (see for example M). If neutrinos were 
sufficiently massive to cluster in our galaxy and make up the local dark matter 
halo (m„ > 20el/), only then would we have / = piocai/m u = 0.4 x 10 7 /m v r e v\- 
This possibility is already rejected by structure formation arguments. With sub- 
eV neutrino masses the enhacement f due to clustering could only be of close to 
1. 

All of the above accelerations are extremely small and beyond the reach of 
any known experimental measurement technology. However, it has been noted 
(by Zeldovich and Khlopov |L8| as well as Smith and Lewin M) that coherent 
scattering of neutrinos from domains of the size of the deBroglie neutrino wave- 
length Xy = 2nTi/p v dramatically increases the scattering cross section. The 
extra factor due to coherence is the number of nuclei in this domain (since the 
nuclear coherence factor is already included), 

N c = —T-P{gr cm- 3 } (K{cm}) • (51) 

For relativistic and for unclustered non-relativistic relic neutrinos with \ u = 
2-Kh/ATy « 0.2 cm this enhancement factor is 

5 x 10 21 ,„. 

A l v{\.$°K} 

For clustered non-relativistic neutrinos with p„ = mv u m 10 m v s e y\ eV 
the wavelength is X v rs (0.12/m„.r e v r }) cm so that the coherence enhancement 
factor is 

10 21 

Nc = A m r ^ P{gr cm ~ 3} ' ^ 

Thus the largest acceleration values with the Nc enhancement factors are 
a R = 2x 10 -34 fp {gr cm - 3} - , (54) 



aNC-N r — 3 x 10 / (m„{ CT }) (T V { lgoK y) P{ gr cm - 3 } 2' (^) 

and finally, for clustered, non-relativistic neutrinos 

aC-NR ~ 10~ fP{grcm- 3 } J' ( 56 ) 

The above discussion applies only to Dirac non-relativistic neutrinos as only 
these have coherent vectorial Z° couplings in the static limit 

A { °] • V^oV- « A%}il> + ip . (57) 
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In the Majorana case we still have coherent ■0 + 75'0^z° couplings. The latter 
are, however, suppressed by j3 v , the ratio of "small" and "large" components 
of the spinor. Recall that for non-relativistic spinors the lower components are 
"smaller" than the upper components by a factor of /3. Without the 75 these 
terms are dwarfed by the leading order terms; however with the 7^75 vertex in 
the static limit these terms are the leading order terms). Thus the analog of 
Eqs. (pS) and (pq) for non-relativistic Majorana neutrinos is suppressed by an 
extra factor of pf ss 10~ 6 . 

3.1 Solar Neutrinos and WIMPS 

Solar neutrinos provide a directional fairly well known source of relativistic 
neutrinos, and it is interesting to estimate their contribution to the accelerations 
we have calculated. The acceleration due to solar neutrinos is 

N av E v 

a = 5>soiar v —r-po(v-A) — ■ (58) 

A C 

For a galium detector, for example, we take $„ ~ 10 n cm~ 2 sec _1 , E v ~ 0.3 
MeV, and ov—gaiium rs 10 -44 cm 2 . Thus for solar pp neutrinos in galium we 
find 

a « 10~ 27 cm/sec 2 . (59) 

This exceeds most previously calculated accelerations, but it is still too small 
to be detectable at present. 

Our main focus here is on neutrino induced forces. However, it is believed 
that much of the cosmological and most of the halo dark matter is made of 
massive weakly interacting particles (WIMPs). 

Let us first estimate the effect of WIMPs if they had the cross section of 
Dirac neutrinos with masses mx ~ 0(100 GeV). Then the nuclear cross section 
with A ~ 100 targets would be large 

o x -a « ^ {m x f A 2 « 10- 3t) m X{GeV} A 2 cm 2 (60) 

and the recoil energies, mxPx/^ ~ O(30 keV) (using mx ~ lOOGeV), de- 
tectable. Indeed such WIMPs have been excluded by direct searches in under- 
ground detectors. Realistic WIMP candidates at the "threshold of detectability" 
have smaller cross-sections by a factor 10~ A . Using the analog of Eq. ( |45| ) 

« Nav 

«WIMP — "PWIMP — J— ^(X-A^xVx 

Nav 2 mi\ 

= n x mx—T-cr(x~A)Vx , (61) 
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and a WIMP density nxmx — Pdark(iocai) — 10 24 gr/cm , we find for A = 100 
and a X -A ~ lO" 30 ^ cm 2 

nwIMpral0 -^,(_^™. m 

Clearly owimp dominates over a very large range of cross sections (with 
A < 10, using Eq. (39) or Eq. (54)) the corresponding accelerations due to the 
scattering of light, locally unclustered neutrinos. 



4 Conclusions 

We have calculated the magnitude of the signals expected for realistic cosmic 
neutrino backgrounds in detectors attempting to measure the mechanical forces 
exerted on macroscopic targets by the elastic scattering of relic neutrinos. We 
examined effects proportional to Gf and G F for both Dirac and Majorana neu- 
trinos either relativistic or non-relativistic. We also estimated the contributions 
to macroscopic accelerations due to solar neutrinos and WIMPs in the galactic 
halo. 

The effect linearly proportional to the weak coupling constant Gf vanishes 
in the case of no lepton asymmetry. With a lepton asymmetry, macroscopic 
accelerations for relativistic Dirac and Majorana neutrinos and clustered non- 
relativistic Dirac neutrinos were found in Eq. (Ell) to be of the order of / • 10 
cm/sec 2 , where f is a density enhancement factor which can be at most about 
17 (with a large lepton asymmetry), and at most one order of magnitude larger 
(i.e. / • 10~ 26 cm/sec 2 ) for non-clustered Dirac or Majorana neutrinos (which 
are most of the relic neutrinos if they are non-relativistic at present). These 
accelerations are at most thirteen orders of magnitude smaller that the smallest 
measureable acceleration of 10~ 12 cm/sec 2 . The acceleration of non-relativistic, 
gravitationally bound Majorana neutrinos vanishes. 

For the effect proportional to G F accelerations of relativistic Dirac and Majo- 
rana neutrinos, taking advantage of coherent scattering effects, were found in Eq. 
( p4|) to be of the order of / • 10~ 34 cm/sec 2 . The accelerations of non-relativistic 
Dirac neutrinos were calculated to be of the order of / • 1Q~ 28 (m v / eV) 2 cm/sec 2 
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in Eq. (55) for non-clustered neutrinos and / • 10 -27 cm/sec 2 in Eq. ( |56| ) for 
clustered neutrinos, while the accelerations of non-relativistic Majorana neu- 
trinos are down by a factor /3 2 ~ 10~ 6 . All accelerations are well beyond the 
smallest measureable acceleration, 10~ 12 cm/sec 2 , mentioned above. 

Additional calculations for the accelerations due to solar neutrinos and WIMPs 
in the galactic halo raise concerns that signals in a detector due to relic neutrinos 
may well be washed out by solar neutrino or WIMP events unless directionality 
can be used to reject them. 
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